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Abstract 

We rewrite the zero-counting formula within the critical strip of the Riemann zeta function 
as a cumulative density distribution; this subsequently allows us to formally derive an integral 
expression for the Li coefficients associated with the Riemann ^-function which, in particular, 
indicate that their positivity criterion is obeyed, whereby entailing the criticality of the non- 
trivial zeros. We conjecture the validity of this and related expressions without the need for 
the Riemann Hypothesis and also offer a physical interpretation of the result and discuss the 
Hilbert-Polya approach. 
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1 Introduction 

The zeros of the Riemann zeta function have been extensively studied in mathematics. Trivial zeros 
of C(s) occur at the negative even integers. In 1859, Riemann [1] conjectured that all non-trivial zeros 
have a real part i. See [2,3] for nice reviews. 

The idea to apply physics to the Riemann zeros is an old one. Hilbert and Polya [4] independently 
proposed that non-trivial zeros of the Riemann zeta function could correspond to eigenvalues of a 
self-adjoint (Hermitian) operator. Subsequent research is suggestive of this spectral interpretation of 
the zeros of the zeta function. The Selberg trace formula [5] encapsulates the duality between the 
lengths of geodesies on a Riemann surface and the spectrum of the Laplacian on the surface. For the 
Laplacian acting on the space of PSL{2, Z) invariant real analytic functions on the upper half plane 
the eigenvalue distribution is related to the Selberg zeta function Z{s) in the same way that Weil's 
explicit formula [6] from analytic number theory is related to the Riemann zeta function. 

Subsequently, Montgomery [7] observed that the pair correlation function of zeros of the zeta 
function along the critical line follows the distribution X{u) ~ 1 — This gives the density of 

spacings of not necessarily consecutive zeros. Dyson realized that the density X(u) specifies the pair 
correlation function of the eigenvalues of a large random Hermitian matrix with a Gaussian measure. 
The statistics of the zeros of the zeta function therefore follow the Gaussian unitary ensemble, a fact 
verified spectacularly through the numerical studies of Odlyzko [8] . 

Over the years the spectral interpretation has been explored by various researchers in an effort to 
bring quantum physics and statistical physics to bear upon this difficult mathematical problem. A 
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large list of references is included in [9] . More recently the connection between physics and the Riemann 
zeros has been explored in the works of Berry and Keating [10] and its subsequent elaborations by 
Conncs [11] and by Sierra and Townscnd [12] regarding the physical interpretation of the formula for 
the number of non-trivial zeros in a given interval. 

The aforementioned formula can be written explicitly in terms of the Riemann ^-function, which 
was also used by Li [13, 14] to establish an equivalent formulation of the Riemann Hypothesis. The 
present article attempts to establish a connection between the physics of the exact formula for the 
number of non-trivial zeros N{T) up to a given height in the critical strip and Li's positivity criterion. 
Our arguments rely on the property of N{T) as a cumulative density function. The intuition and our 
discussion find their wellspring in theoretical physics. 

The note is organized as follows. First, we review Li's Criterion, with an emphasis on requisite 
properties of the Riemann ^-function. In the central section of the paper, we show how using N{T) 
and its derivative as eigenvalue distributions in the critical strip would imply the positivity of Li's 
coefficients. Finally, we interpret the result from a physicist's perspective. 

Nomenclature 

Throughout this paper, we will adhere to the following notation. 



^{z) The Riemann ^-bmction, £,{s) := ^s{s - l)7r-ir(s/2)C(s). 

S The critical strip, 5 := {z | < Rc{z) < 1}. 

iS+ The upper critical strip, S+ := {z \ < Re{z) < 1, Im(z) > 0}. 

jC The critical line, jC := {z \ Re{z) = ^}. 

jC+ The upper critical line, C := {z \ Re{z) = 5, Im(z) > 0}. 

N{T) The number of zeros of ({z) or ^(z) inside <S+ up to height lm{z) < T. 



Pj,k = (^j,k + il^j A zero of ({z) or ^{z) in S+ indexed by j, k G Z>o; Hj+i > Hj, (yj,k+i > cTj.fe- 

2 Li's Positivity Criterion 

Let us, perhaps more for the sake of notation, first remind the reader of the statement of the Riemann 
Hypothesis [1,15-17]: 

The analytic continuation, from z G ]R>i to the whole complex plane z G C of ({z) := 

00 

J2 has all its zeros in the critical strip <S := | < Re(^) < 1} lying on the critical 

71=1 

line jC := {z \ Re{z) = 5}. 



3 



Now, C(z) obeys a remarkable functional equation C(z) = 2^7r^ ^ sin(i7rz)r(l — z)C(l — z), which 
inspired Riemann to define the ^-function, which will be central to our discussions: 

e(z):=iz(z-lK-tr(|)C(z) . (2.1) 

There are many advantages to considering ^(2) instead of C{z)i which we now summarize. 



2.1 The Riemann ,^-Function 

First, note that ({z) has trivial zeros at all negative even integers, which are conveniently cancelled by 
the corresponding simple poles of the F-function. Hence, ^(z) has only non-trivial zeros, all located 
within the critical strip, and affords an elegant Weierstrafi product expansion, known as the Hadamard 
product: 

CW=e(0)n(l--)' (2.2) 

pes ^ 

where ^(0) — ^, and p are the non-trivial zeros of C,{z) in S since none of the other factors in the 
definition of ^(z) vanishes therein. It is an obvious but important fact that S extends to both the 
upper and the lower half-planes and that all the p occur in conjugate pairs above and below the real 
line. We shall denote S+ as the critical strip above the real line and correspondingly as the upper 
critical line. 

Second, the functional equation becomes particularly simple: 

az) = a^-z). (2.3) 

Indeed, the symmetry about the critical line of Re(z) = i becomes manifest. Hence, not only are the 
zeros symmetric about the real axis, they are also symmetric about the critical line. The Riemann 
Hypothesis postulates that all the zeros in fact lie on £. 

Third, consider the conformal mapping z <—> ^^j, which maps to the boundary circle |z| = 1 
of the unit disk and the entirety of the left half of 5+, i.e., {z | < Re(z) < i}, to the interior of the 
open unit disk. Consider, therefore, the function 

0(z):=e(^)-e(-^), (2.4) 
z — 1 1 — z 

where we have used (j2.3p in the second equality. Since all zeros are symmetric around the critical line, 
to consider the zeros of within S therefore amounts to considering the zeros of (/)(z) within the open 
unit disk. Indeed, the Riemann Hypothesis would require that there be no such zeros and indeed that 
all the critical zeros are on the boundary circle. Since the number of zeros of an analytic function 
/(z) in a region R is equal to J^^ dz we have that 

PROPOSITION 1 The Riemann Hypothesis is equivalent to being analytic (i.e., holomorphic 

and without poles) within the open unit disk \z\ < 1. 
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Fourth, a beautiful exact result relates the number of zeros up to height T within S, and the 
^-function: 

N{T) := #{p e 5 , < Im(p) < T| ^{p) = Q} = - Imloge(i + iT) . (2.5) 
Briefly, this can be seen [15] as follows. We have that, since there are no poles in the critical strip, 

^(-^) ~ 2^1 IdR'^'^ J(sY ^ In^^ Isr'^'^ (logC(s))'i where R is the rectangular region in the critical 
strip up to height T, strictly, R = {z \ —e< Re(z) < 1 + e, < lm{z) < T}. (Note that no zero passes 
through the hne Im(z) = T.) By the symmetry of we have that N{T) = 2^ Im ds (logf (s))', 

where C is the right half segment of dR from i to i + iT, and whence (|2.5p since log^(s) is real at i. 

Extensive work has been done in estimating (|2.5p . We have, using the definition (|2.ip for the first 
equality, that 

NiT) = -^{T) + 1 + - Imlog C(l + zT) ~ ^ log ^ - ^ - I + 0(log T) , (2.6) 
TT TT 2 27r 27r 27r 8 

where historically ^{T) = Imlogr(|r+ j) — -j logTr is known as the average part and ^ ImlogC(^ + 

iT) C'(logr) is the fluctuating part around the essentially TlogT growth of the former. This is the 

inspiration behind classical (the average) and quantum mechanical (the fluctuation) interpretations 

of the critical zeros [4, 10]. 

We note that N{T) is a real step functiorl*], increasing by unity each time a new critical zero is 
encountered: 

oo 

NiT)^ J2O{T-lmp) = Y,aJ0{T~^^,) . (2.7) 

It is vital to explain the above in detail. We can explicitly write the upper critical zeros as pj,k — 
<^j,k + il^j, indexed by j,k e Z>o, where aj,k G (0, 1) and fij G M>o. The zeros are ordered so that 
Pj+i > Pj. Crucially, as we do not assume the Riemann Hypothesis, the real part of p need not be 
i. Furthermore, we do not assume simplicity of the zeros, which is as yet also unknown [18]. If, for 
example, we have a double root, we explicitly count this twice. By the functional identity, zeros of 
the ^-function not on the critical line are paired within the upper critical strip: if cr + i/i is a zero, 
then so is (1 — cr) + ip. The index k = I, . . . ,aj enumerates the zeros with the same imaginary part. 
We order the zeros so that <7j.k+i > '^j.fc- The aj then counts the number of zeros with imaginary 
part pj including the multiplicities of the roots. If, for example, there are a pair of simple roots 
with imaginary part pj, then in an cpsilon interval around pj, the counting function N{T) jumps by 
two. Contrariwise, if p G is the only simple root with imaginary part pj, then N{T) jumps by 
one. We know that the number of roots with imaginary part in the interval (0, T) is finite. Indeed, 
the asymptotics of the expression are given in (|2.6p . In summary, N{T) is a strictly increasing step 
function as we move up in height T regardless of the Riemann Hypothesis. 

f 1, x>a 

* We use the standard notation that the step function is defined as 9(x — a) = < . Its derivative is the 

[ 0, X < a 

delta-function 5{x — a) = 9'{x — a). Importantly, the integral of 5{x — a) over any finite interval containing a on the 
real axis is normalized to equal to unity. 
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The first few numerical values for /ij=i,2,3,4,... are approximately 14.134725142, 21.022039639, 
25.010857580, 30.424876126, .... It is known that the first 10" zeros lie on the critical line [19]. The 
Riemann Hypothesis is, of course, the statement that aj = ^ for all zeros. 

The rewriting of N{T) lends itself to a wonderful interpretation: N{T) is a cumulative density 
function defined over the critical strip. In other words, its derivative, p(T) :— N'(T) is a density func- 
tion of distributions. That is, one could conceive of a physical system whose energy levels (eigenvalues 
of the Hamiltonian) are thus distributed; this is along the school of thought of Hilbert-Polya [4]. Of 
course, the resulting eigenvalue density is highly non-smooth, but is, rather, a sum of delta- functions: 



(2.8) 



2.2 Li's Criterion 

Another striking property of the ^-function was noted by Li [13] not too long ago. Let {fc„}, for 
positive integers n, be defined by 

1 d" 



kji : — 



-,"-1 



{n - 1)! dz" 



(2.9) 



then we have that 



THEOREM 2.1 [Li's Criterion] The Riemann Hypothesis is equivalent to the condition that fcn > 
for all n G Z>o. 



Li also showed two equivalent ways of writing these numbers, namely 



E 

pes 



where we recall that 



(2.10) 
(2.11) 

(2.12) 



Therefore, (|2.10p means that fc„ are simply the Taylor expansion coefficients of (log^(Y3^))'. We 
stress that in order for this Taylor expansion to make sense, the radius of convergence must be such 
that there are no poles within the enclosed disk (within the unit disk) thus defined. That is, we cannot 
allow a zero to occur for the denominator. Let a zero he p := cr + i p, then the first zo value for which 
~ p will restrict the radius of convergence. Whence |zo| = [13], to be safe, a worst 

bound is assumed, where a is taken to be 1 and |zo| = j- 

On the other hand, (|2.1ip is also important [13, 14]; note that the sum is real because all critical 
zeros p occur in conjugate pairs. Furthermore, note that the modulus 1 — ^ can neither exceed nor 
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be exceeded by unity since the summand 1 — [1 — - ] in the former case grows polynomially in n 



and the sum diverges, while it tends to 1 in the latter case and the sum also diverges. Indeed, if the 



Riemann Hypothesis held, then for the j-th pair of conjugate critical zeros 



1 



2±'lfJ.j 



which is exactly equal to 1 and the sum in (12.111) converges. Moreover, the unimodularity allows us 



to set 1 — T-;- — 



exp(±i0j), whence we would have 



kn = ^(1 - exp(m6'j)) + (1 - exp(-m6'j)) = 2^(1 - cos{n9j)) > 
Now, for the converse. Li explicitly calculated the coefficients a, in (|2.12p to be 



(2.13) 



a, =4^(^^._^)i ^°°dx [a;3/V(x)]'Qlogxy[l + (-l)%-i/2]eR>o , (2.14) 



where 'ip{x) := J2 ^ 



Using (|2.10p and (|2.12|) , one has the recursion relation that 

n-l 

kn — ^ ^ ^ kj dji^j • 

J = l 



(2.15) 



If all kn > 0, then the above recursion would imply that fc„ < nan for all n G Z>o, whence for z in 
the unit disk, 



< |fc„z" i| < ^na„ \z\' 



0'(kl) <oo 



(2.16) 



That is, is analytic in the open unit disk and by Proposition[Tl Riemann Hypothesis holds. Thus, 
the necessary and sufficiency together implies Theorem 12.11 

Explicitly, one can Taylor expand to find the Li coefficients. The first few are: 



fci 



1 



_(2 + 7 + log-) , 



, 7r2 1 
k2 = 1 + 7 - 7 + ^ + log — - 271 , 
8 47r 



(2.17) 



1 - 37^ + 7"^ + — + - log 



1 1 

4^ ^ 16' 



V^2(l)-67i+7 - + 371 +-72 -C(3) 



where 7 is the Euler-Maschcroni constant, 7^ are the Stieltjes constant Jl|, and tj^niz) — gprrr logr(z) 
is the polygamma function. Indeed, we present some of the first numerical values of these coefficients 
in Table □ 



t We recall the definition of the Stieltjes constants. Consider the series expansion of 1^(2:) about z = 1: 

1 °° C— 11" 

(logfc)" (logm)"+i' 



7n 



lim 

m — 'oo 



n + 1 



The 7„ are then the Stieltjes constants. We can as well define 

d" , ^ {-!)" n! 

7n = -r-rC(^) 



dz" 

In particular, 70 = 7, the Eulor— Mascheroni constant. 



3 Density Function for the Distribution of Critical Zeros 



Thus armed with all necessary ingredients, let us proceed to rewrite the Li coefhcients in a suggestive 
form. We will do so formally in two ways and point out the subtleties involved at the end of the 
manipulations. 



3.1 A Logarithmic Expansion 

First, integrating (|2.10p and using the definition as well as the functional equation (|2.4p . we have that, 
for z in the unit disk, 

-.00-. 00 , 

loge(^)-loge(T— / dzz"-i = -log2 + ^^z" , (3.1) 

n=l n=l 

where — log 2 is easily checked to be the constant of integration. 

We note, however, that the above expansion has radius of convergence strictly less than unity, 
whereby excluding the unit circle, to which, crucially, C is mapped under our conformal transformation 
z I— > Therefore, it is imperative to analytically continue. Let us start with p.ip and rewrite the 
expansion about the point z = — 1. Note that the region of convergence for this is an open of circle of 
radius two centered about z — —1, which, in particular, encloses the entirety of the closed unit disk, 
except the point z = 1, where there is a pole for (p{z) = ^{j^)- Therefore, we have that 

log^(Y^) = log- + ^^z" = 6o + ^fe„(z + l)" , (3.2) 

n—l n—1 

where we have written the two expansions together for comparison. Indeed, as with the Li coefficients, 
all bn are clearly real. Note that 

_rfi)(:fi) 1 / 1 1 \ 

bo:=log g ' . 61:= g (^-log^ + V^o(^) + 2[logC(-)]'j =0, (3.3) 

as well as the curious identity 62 = &3- 

Expanding (z + 1)", we readily obtain an expression for the Li coefhcients in terms of the new 
expansion coefficients 6„: 

00 ^ . s 

J 



We now wish to solve for the kn explicitly and demonstrate positivity. 

We first see, using the counting formula (|2.5[) and the expansion (|3.2p . that 



iV(/,) = i Im log C(i+z/.) = f:^Imf 1^ + 1)" , (3.5) 

^ " ^-^ TT \2fJ. — I J 



n=l 
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where we have used the substitution ^ in = -r—i or z — Then, since 

Im ( + iV = , ^y",,. sin(ntan-i — ) = 2" cos" 9 sm{ne) , cosO , , 

(3.6) 

we have that 

oo 

TrN{n) = J2 ^"2" cos" e sm{n0) . (3.7) 



Next, we apply the integral identitjcl 

7r f Ti\ 

d6l cos"6lsin(n6l)sin(2m6l) = — — , m, n e Z>o . (3.8) 

2"+^ \m J 

Therefore, multiplymg both sides of (13. 7p by sin(2TO6') d^ and integrating, we find that 

/ ' de* TiNU) sin(2TO6l) = V 6„ ( " ) - . (3.9) 

Whence, 

y =4 d0iV(Ai) ;7™_i(cos(20)) sin(20) , (3.10) 

where we have used the definition of the Chebyshev polynomial of the second kind, that 

cos :^ : 7^ . 3.11 
sin0 

Finally, using that cos(26') = cos^ 6 — sit? 9 = f^^^TT' ^^^^ sin(20) d9 = — 2cos0 d(cos0), and that as 
fi proceeds from to cx), subtends an angle from -I to 0, we find: 

gbj") =8 rdM^(/.) f/_,(l4^)^^^-^ . (3.12) 

^1 VW Jo V + 1 ^4,^2 + 1(4/2^ + 1)2 

However, recognizing the left hand side from (j3.4[) . we arrive at the conclusior0 that for all m e Z>o: 

r°° u 4u2 - 1 

k^^32rnj^ ^^-J^ N{^,) U^^,{J^^) . (3.13) 



* We can show this using the elementary identity that 2" cos" 9 sin(nf?) = (") sin(2p6'), as well as the orthogonality 

condition that d9 sm(2p9) sin{2m9) = —jnSpm, p,m £X. To our knowledge, 1 13. 8| I first appeared in D. Bierens de 
Hahn's Nouvelles tables d'integrales definies, Leiden: P. Engels (1867), published one year after Riemann's death. 

t Curiously, if we did not worry about the region of convergence and proceeded formally with 1)3. we can obtain 
the same result. Indeed, we can rewrite the function 112.51 as, for fi G IR>o> 

iV(M) = i Im log C(i + i/.) = £ ^ Im r + ' ' 



^2fj, — i 

where we have set i + i/^ = or z = for the expansion and where we have used the fact that fc„ are all real. 

Clearly, ^^^^ is unimodular and resides on the unit circle. Using l|3.11ll . we can simplify the expansion parameter: 

\m(— =Ime ^I^T^' = sin tan~^ ra = j7„-i(-h- . 

Next, we recall that the Chebyshev polynomials prescribe an orthonormal basis over [—1,1]: 

/ Ax Un{x)Um(3:)\/ 1 - = —Smn , 
J-1 2 
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This is beautiful and remarkable integral expression for the Li coefficients. We numerically evaluate 
some of the first of these numbers and find them to be indeed close to the required values of Li. This 
is reassuring, and the results are presented in Table [TJ 

Let us now proceed to simplify the integral (|3.13p . First, recall that Chebyshev polynomials of the 
second kind are related to the Chebyshev polynomials of the first kind by the equation 

dx Unix) = ^— T„+i(x) . (3.14) 
n + 1 

Next, integration by parts J^^ dn f (fJ.) g' (tJ,) = [fifJ-] 9i^J■)]a ~ d/i f'Ui)g{^) suggestively compels us 
to define 

/(^) := N{^,) ■ g'i^) := ^^g^^^, t/„._i(g^) =^ g{^,) = 2(r™(g^) - 1) . (3.15) 

Note that we have added —2 as the integration constant for g(/i) for reasons shortly to be clarified. 
The integral thus becomef^ 

2 1 -\°° poo , 2 



27V(m)(T„,(J|^)-1) 



- 2 / d^ p(^)(r,„(^JJ^) - 1) , (3.16) 



J 

where, importantly, we have used (|2.8p for the explicit form of p(/i) and the defining property of the 
delta-function. Next, note that A^(0) = since the first critical zero does not occur until /i > 14, and 
that for any fixed m. 



hm iv(^)(r™(j4-T) - 1) ~ ™'(rTTT) = ^ ■ (3-1^) 

Thus, the boundary terms vanish, and we have that 



fc„ = 2^aj (l-r™(-4— )) , meZ>o. (3.18) 

Recalling the definition of the Chebyshev polynomial of the first kind, that Tm(cos6') :— cos{m6) 
and remembering that all aj, regardless of actual value, are strictly positive integers, each summand 



which, when transformed to x := "^^-i , ] , becomes 



Jo ^ ^"^i;?TT^^'"^5;?TT^4M2 + 1(4^2 + 1)2 " i"^™" • 

Therefore, this allows for the inversion of the formula for N{ii) above by integration of both sides, and we arrive once 
more at 113.131 1. 

One could also perform the integral without recourse to delta-functions and generalized analysis. The expression 
ll3?T3l l becomes, by ifZTj l. 

/'°° fj, ip? — 1 ^^i, r°° 16/j 4/^2 _ 1 

= 32mX:a,y^ " + 1)2 ^ VTT^ = g 7,^ ''^ (4^2 + 1)2 "^-^^ ^TT^ 

1 oo 4/(2 



' — ^ ]^ -■- ^ — 4/j, — 1 



-T- ,-1 
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in (I3.18P is non-negative. In fact, because | < 1 strictly for any imaginary part as we know the 
first is at least 14, 

4^2 _i 

Tm{-r-^ 1") ^ given rn e Z>o and for all j G Z>o , (3.19) 

4:flj + 1 

so each summand is actually strictly positive. In conclusion, this would imply that all Li's coefficients 
km > 0. Indeed, we have numerically evaluated the first few km coefficients given in (IS.lSp . using the 
first 10'* critical zeros, and see that they are very close to the results obtained from (|2.10p and (|3.13p . 
and moreover constitute an increasing sequence of strictly positive numbers. In particular, the above 
discussions would force all critical zeros to lie entirely on £. In this case, treating any multiple roots 
distinctly, we see that ()3.18p reduces precisely to (|2.13p . which saturates the lower bound. 

We seem to have arrived the positivity of the Li's coefficients without the assumption of the 
Riemann Hypothesis! This, of course, is not quite true: the criticality of the zeros of the ^-function 
has in fact implicitly been invoked^] Let us re-examine the expansion p.ip . If a zero of f were 
ever encountered, a branch cut must be carefully chosen because of the logarithm. Now, in setting 
+ ill = in ()3.5p . we have used the fact that /i is real, whereby making z on the boundary of 
the unit disk to which the critical line is mapped. That is, z — is always unimodular for real fj,. 
Therefore, only when no zeros at all lie in the interior of the unit disk, to which the region inside S 
but not on C is mapped, do we not need to carefully find the branch. This, then, would require the 
Riemann Hypothesis. Nevertheless, we have formally arrived the integral formula (|3.13p which does 
not seem to depend on the precise location of the zeros and the subsequent choice of branch cut. We 
will discuss possible implications of this later. 



3.2 A Delta-Function Expansion 



Can we evade the problem of logarithmic branch cuts? Let us return to the original definition of the 
Li's coefficients from (|2.10l) : 



' ^ ' ^^l-z' n=0 n=0 

Thus defined, the relation between the k and b coefficients in (13.41) still holds. 



(3.20) 



Now, let us, with the aid of some non-standard analysis, reinvestigate the counting function N{ii) = 

i Im log ^ ( i + i /i) =^Oij 0{iJ, — fij), where /i j is the imaginary part of the j-th zero in 5+ , regardless 

j 

of the Riemann Hypothesis, and aj counts potential multiplicities. Upon taking the derivative, this 
gives us a sum over delta-functions, or. 



Sill- fij] 



■Re 



C(i+*M) 



(3.21) 



We are grateful to R. Heath-Brown and H. Bui for pointing this out to us. 
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n 


kn by expansion 


kn by integral 


DifTcrence 


kn by sum 


Difference 


1 


0.0230957 


0.0235290 


-1.876% 


0.0229610 


0.5832% 


2 


0.0923457 


0.0940796 


-1.878% 


0.0918069 


0.5835% 


3 


0.207639 


0.211543 


-1.880% 


0.206427 


0.5838% 


4 


0.368790 


0.375740 


-1.884% 


0.366635 


0.5844% 


5 


0.575543 


0.586417 


-1.889% 


0.572175 


0.5851% 


6 


0.827566 


0.843252 


-1.895% 


0.822717 


0.5860% 


7 


1.12446 


1.14585 


-1.903% 


1.11786 


0.5870% 


8 


1.46576 


1.49376 


-1.911% 


1.45713 


0.5881% 


9 


1.85092 


1.88645 


-1.920% 


1.84001 


0.5895% 


10 


2.27934 


2.32332 


-1.930% 


2.26587 


0.5910% 


11 


2.75036 


2.80373 


-1.940% 


2.73406 


0.5926% 


12 


3.26326 


3.32695 


-1.952% 


3.24386 


0.5944% 


13 


3.81724 


3.89222 


-1.964% 


3.79448 


0.5964% 


14 


4.41148 


4.49869 


-1.977% 


4.38508 


0.5985% 


15 


5.04508 


5.14550 


-1.990% 


5.01477 


0.6007% 


16 


5.71711 


5.83170 


-2.004% 


5.68263 


0.6032% 


17 


6.42658 


6.55632 


-2.019% 


6.38765 


0.6057% 


18 


7.17248 


7.31832 


-2.033% 


7.12884 


0.6085% 


19 


7.95374 


8.11664 


-2.048% 


7.90512 


0.6114% 


20 


8.76928 


8.95012 


-2.062% 


8.71540 


0.6144% 



Table 1: The Li coefficients evaluated numerically. Here, "by expansion" means the Taylor expansion in 
(|2.10p , "by integral" means the integral identity in (|3.13p . with cut-off 10* for the upper bound, and "by 
sum" means the summation over Chebyshev polynomials in (|3.18p . with cut-off as the first 10^ critical 
zeros. We have also juxtaposed the percentage errors with respect to the expansion values for comparative 
purposes. 

where we have exchanged the derivative with the Im. 

Again, using the conformal transformation i + i/i = or z = we have 
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Now, we must substitute the following relation: 

= (2" cos" 61) (4 sin^ 9) cos{n0 + V) (3.23) 
= 2"+2 cos" e sin^ e cos((n + 2)61 - tt) , 

where we have used, upon defining the angles, 



1 4/i 

tan0 = — , tanV- = — tt =^ 29 ^ tp + tt . (3.24) 

2fi Aji^ — 1 



Simplifying cos((n + 2)6 — tt) = — cos((7i + 2)6'), wc find that 

oc ^ oo 

^ Uj 6{fi- fi-i) = '^{n + 1) bn+i 2"+2 cos" sin^ 6 cos((n + 2)61) . (3.25) 



TT 

j=l n=0 



For each positive integer m, we shall hit both sides of p.25p with Asm^{m6{ii)) dii and integrate 
the variable fi on the half-line [0, oo). Recall that 2sin^(m^?) = (1 — cos{2m9)). Thus, the left hand 
side becomes 

poo °° °° 

/ d/x^aj 5(^-Mj) 2(l-cos(2m6l(^))) = 2 ^ (1 - T„,(cos(26'(^j))) 

°° 4u2 - 1 

= 2^a, (l-T„,(-^4— )) , (3.26) 

j=i ^ 

where we have used the defining relation of the Chebyshev polynomial of the first kind. 

Now, we do the same on the right hand side. From (|3.24[) . ^ e [0,cx)) implies that 6 G [§,0). It is 
easy to check that d/i = — ^^.^^^j g d6>. We as well note the definite integral identities 

/ ' d6i cos"6l cos((7i + 2)6') = 0, n e Z>o , (3.27) 
Jq 

{-JL-( ) n> m-2 
2" + - Vm-i; ' ' (3 2g) 

, 71 < m - 2 . 

Putting these pieces together, we find 

V (n + 1) bn+i 2"+2 / d6' cos" 9 cos((n + 2)61) (1 - cos(2m6l)) 



, m — 1/ \m 

n—m—l ^ ^ n—m 



J2 bn+iin + l)i " ) =m ^ 6„ P ) = fc„ . (3.29) 



Hence, we recover p. 181) 

oo 

_ = 2 Vn-. n - T„,l 



°° 4u2 - 1 

fc™ = 2^aj (1 - Trni^J^-^)) , m e Z>o . (3.30) 
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Again, we seem to have arrived at an expression which imphes Li's positivity without assuming 
Riemann Hypothesis. However, re-examining p.20p reveals a subtlety. The radius of convergence is, 
as pointed out in [13], for \z\ < j. However, we have equated z = j0j P-^^p . because /i is real, and 
\z\ = 1, whereby again making the expansion and subsequent derivations only formal. Of course, one 
could argue that the delta-function involved is a generalized function and should be neither confined 
to nor validated by convergence, since it is, after all, itself a divergent quantity. Nevertheless, one 
needs to rigorously define and manipulate the divergences carefully. 

3.3 Remarks 

Here we wish to collect some comments regarding our main formulas for the Li coefficients: 

/■oo ^ 2 _ °o 4„2 _ ^ 

K^32nJ^ d^^ ^^^/^ N{^^) ^ 2^ a, (1 - ■ (3-31) 

j — 1 J 

These formulas have been derived by formal manipulations of divergent series and then have been 
checked numerically. The crucial question is whether they are true in a rigorous mathematical sense. 
In particular, there might exist a contour deformation that takes the integral formula to the convergent 
region and leaves it unchanged. We leave that question as open, and list these formulas as hypothetical 
assertions to be proven rigorously. Li also shows that the Riemann Hypothesis for the Dedekind zeta 
function C-R"!*) is as well equivalent to the non- negativity of a sequence of real numbers [13]. We 
therefore conjecture that similar formulas to the expressions for fc„ apply for the coeSicients k„ in the 
series expansion of (log^^l jttj))'- 

Another question is whether these formulas assume the Riemann Hypothesis in some subtle way. 
The derivations we have presented indeed suffer from this affliction, but it may be that an alternate 
derivation of this expression sidesteps the issue. We will offer a few further comments on this score. 
The crucial observation is that both the counting function N{fi) and the argument of the Chebyshev 
polynomials that appear in the integral formula for the Li coefficients know directly about factor of 
i only from the functional equation for the Riemann ^-function. One might also question whether 
the unimodularity of the trigonometric representation of the Chebyshev polynomials is the signal of a 
secret assumption of the Riemann Hypothesis, but once again we note that this unimodularity relies 
on the argument ^ + i/^ of the counting function N{ii) determined by the functional equation for the 
^-function. This tempts us to think that the final integral and summation formulas for fc„ only know 
about the functional equation for the Riemann function ^. But this remains to be rigorously proven. 

Finally, we wish to offer a couple of comments about the possible geometric and topological meaning 
of the integral and the summation formulas for Here we wish to compare the integral and the 
summation formulas for the Li coefficients to the remarkably deep structures uncovered in the profound 
proofs of the Riemann Hypothesis for the case of finite fields by Weil, Deligne and others [20]. (An 
insightful simrmary of Weil's work is presented in the recent book [11], especially in chapter 4.) The 
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work of Li is also related to the seminal work of Weil by Bombieri and Lagarias [14] . We wish to note 
that the integral and summation formulas for fc„ discussed in this paper capture, at least heuristically, 
some of the crucial aspects discussed in these papers and books. 

In particular, the equality of the integral and the summation formula for the Li coefhcients reminds 
one of Weil's explicit formula [6]. This explicit formula has an interpretation in terms of an index 
theorem with deep meaning in algebraic geometry. The integral formula for the Li coefficients could be 
understood as an index formula of the Atiyah Singer type, which might have a profound topological 
interpretation in the case of the distribution of the Riemann zeros as discussed in Connes and Mar- 
colli [11]. For example, the counting function appearing in the integral formula for A;„ might be seen 
to correspond to some generalized Todd genus and the Chebyshev polynomial to some generalized 
Chern character. The summation formula for the Li coefficients on the other hand can be recast as an 
alternative sum, which is very similar to the Lefschetz fixed point formula which figures prominently 
in Weil's seminal work and in more recent topological discussions of the Riemann Hypothesis, as once 
again reviewed in [11]. Finally, even the argument of the Chebyshev polynomials that appears in the 
integral and the summation formulas for fc„ is suggestive of the existence of some non-trivial vector 
bundle structure over the space of integers, where the factor t^i,] reminds one of the value of the 
Higgs field vortex configuration which determines the self-dual connection on that vector bundle. (See 
the Section 7.6 of Chapter 4 of Connes and MarcoUi [11] for comparison.) The summation formula 
then appears as the formula for the action of the compact (topological) Abelian Ifiggs theory, as 
prominently discussed in the case of compact quantum electrodynamics, as prominently discussed in 
the case of compact quantum electrodynamics [21]. 

We intend to return to these fascinating structures in future work. 

4 Physical Perspectives 

As we have noted in the Introduction, attempts to use a physical argument to explain the Riemann 
Hypothesis date back to Hilbert and Polya [4]. Were one to find a matrix 

H=yI + i.T, (4.1) 

whose eigenvalues are all the zeros of ^(z) and if in addition the matrix T were Hermitian. then all the 
zeros would lie on the critical line C. Therefore, one would wish to find the appropriate operator H 
and its corresponding eigenfunctions, or, in physical terms, the Hamiltonian and its wave functions. 

Berry and Keating [10] pioneered a method to interpret N{T) as a density of states. However, 
they focused primarily on the "average part" of N{T), viz., the part that exhibits TlogT growth. 
We seek to use Li's Criterion in concert with our approach to further this interpretation. The full 
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expression N{T) = i Inilog^(i + iT) inspires us to regard the ^-function as a wave function, 

= r(z)e'^'5(") , (4.2) 

where S{z) is the standard quantum mechanical phase and r{z) is the modulus. Subsequently S{z) = 
;^Imlog^(z) indeed becomes N{T) for z = ^ + iT. In the Hamilton-Jacobi formalism of classical 
mechanics and its standard quantization, S{z) in (|4.2p is simply the classical action. We know that 
N{T) is an integer. Therefore, the action S must be quantized, i.e., S = hn for ri G Z, and S is 
therefore integral in units of h. This requirement is, of course, the Bohr-Sommerfcld quantization 
condition. 

Furthermore, using (|3.5p and (|3.14p . as well as (|3.18p . we have that, for the classical action (i.e., 
dropping the factor of h), 

1 °° 4z2 - 1 oo / 4„2 _ ]^ \ 

Siz)^-J2kndyn , 2/„:=r„(-^— ) , fc„ = 2^ a, 1 - T„(-^^) . (4.3) 

This way of writing an action is well known from Hamilton-Jacobi's action-angle theory. It means we 
have a classical system of an infinite but countable set of particles indexed by n, whose generalized 
coordinates are y„ and conjugate momenta are fc„, with an energy z. Therefore, we propose that the 
Li coefficients should be regarded as generalized momenta and are the action variables. These must 
be positive in order to have a stable classical system. If they were negative, the system would be 
unstable. If they were complex, the system would be metastable. Hence, the Li Criterion implies that 
our dynamical system is classically stable. 

What can one say about the quantum mechanical picture? Let us investigate the secular equation 
of the operator in (|4.ip . It should by construction have roots corresponding to the zeros of (,{z), that 
is, 

f(z) =det(z-I-i/) . (4.4) 

To arrive at the appropriate matrix description, we start by considering the expansion (|3.ip . which, 
combined with (|3.18p . gives 



1 — s n 



n—l j—1 L ' 



J 



(4.5) 



Put cosOj :— ^Jij^i ■ Using the definition of the Chebyshev polynomial and reversing the sums, we 
obtain 

1 oo oo . . 

log^(-— ) = -log2 + ^2a,^ cosine,)) (4.6) 



i=i 



log 2 + ^ aj [-2 log(l - s) + log(l - s e^^^ ) + log(l - s e 
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Substituting z = and using that e^*^J = 1 — (i ± i/i^) ^, we find 

which is simply a refined version of tlic Hadamard product, a reassuring checkj*] Note that each factor 
in the product obeys the functional identity. 

Equation (14. 7p is to be compared with (|4.4p . We see that there exists a matrix such that its 
eigenvalues coincide with /Xj, which is consistent with the Hilbert-Polya picture. As the expression 
vanishes for zeros on the critical line, the refined Hadamard product is already in the form of the 
determinant of an infinite matrix, and the eigenvalues of T are real. Note the appearance of the i 
prefactor in (|4.7p . In a quantum mechanical system, it would be tempting to interpret this as the 
vacuum energy term, but we will give a different interpretation in what follows. 

Knowing the energy levels (eigenvalues) of the Hamiltonian H does not allow us to pin down what 
the latter actually is; the previous discussion simply shows that it must exist. We do know, however, 
that the partition function is 

00 

Z{ii)^lr{e-0")=Y,e-^^^ , (4.8) 

for the inverse temperature /3. This object would define the statistical mechanics of Ricmann zeros. 

What can we say about the actual Hamiltonian? Our discussion here is necessarily heuristic. The 
physical interpretation of the Li formula suggests a system in which the momenta (p) and coordi- 
nates (g) are linearly related. In the simplest case, this would, assuming factorization, correspond 
to the Hamiltonian for an inverted harmonic oscillator H — — 9^), which, for example, features 
prominently in the non-perturbative definition of the two-dimensional critical string theory in terms 
of Matrix models. Upon a further canonical transformation such a Hamiltonian would amount to 
H^QP. 

Berry and Keating [10], Connes [11], and Sierra and Townscnd [12] among others, have looked at 
precisely such a Hermitian quantum mechanical Hamiltonian. For convenience, let us work in units 
where = 1 . We have 

H^i(gP + PQ), P--^. (4.9) 
The Q eigenfunctions of this scaling operator in the Q, P phase space are: 

H i;E{q) = E ijE{q) , ^siq) = (4.10) 

q2 



The product formula l|4.7| l can also be obtained by using the integral formula l|3.13ll derived in section 3. Thus, 



We exchange the order of the integral and the sum and repeat the steps from the previous computation to obtain (14.71 . 
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with some appropriate normalization constant A. The momentum state wave function is the Fourier 
transform: 

^Eip) = r dq e-^^-y.Eiq) = tA^^'"'^^!^^^, ■ (4.11) 

A few further comments are in order. 

• First note that the ^ appears from the symmetrization ^{QP + PQ) in the Hamiltonian. It is 
the same ^ that appears in formula (j2.5p for the number of zeros in the upper critical strip with 
imaginary part less than T and in (|4.7p . We see that this factor is different from the vacuum 
energy, even though it does conspire with the canonical commutation relations (the reason for 
the vacuum energy) to produce the ^ + iE combination. 

• Classically QP is the simplest Hamiltonian that yields the generator of the Mellin transform, 
i.e., a power, or scaling function, and ^ ~ iE appears as the scaling dimension. That is to 
say, the \{QP + PQ) Hamiltonian is selected as the required quantum scaling operator with 
eigenfunctions given by generators of the Mellin transform. 

• The importance of the Mellin transform is indicated in the deep study of Bombieri and La- 
garias [14] regarding Li's positivity condition. In particular Bombieri and Lagarias show that 
Li's Criterion is simply the consequence of the famous Weil Criterion [6], which explicitly re- 
lies on properties of the inverse Mellin transform, namely the fact that the convolution of the 
inverse Mellin transform of a function and its complement (that is to say, its "dual" under the 
s I— > 1 — s map) is equal to the sum of the inverse Mellin transforms. The QP Hamiltonian (even 
classically) enjoys this feature, provided that we identify in the formula for the Li coefficients 
Q„ — e™^J and P„ = e~*"^^. The Hermitian \{QP + PQ) Hamiltonian is the quantum coun- 
terpart of the classical scaling Hamiltonian, which leads to the power law eigenfunctions that 
generate the Mellin transform. 

We observe that there is a beautiful identity, valid for all Re(z) > 0, and in particular, within the 
critical strip, that C,{z)T{z){l — 2^^^) = /q°° da; frqrf- Together with the duplication formula for the 
F-function, ^/^^(z) = 2^"ir(f )F(f + i), one can write 

fca z-1 °o ^cx; 7(7— A\ 

Jo e- + l Jo ^— (2-_2)r(^) 

(4.12) 

This tells us that the ^-function, up to a factor which does not vanish in the critical strip, is the Mellin 
transform of a Fermi-Dirac distribution. Note that, in the second equality we have trivially inserted 
the infinite geometric sum of (— l)"e"^ in order to compare with (|4.1ip . 

Thus, even though the momentum space wave functions of the ■^{QP + PQ) scaling Hamiltonian 
contain the F-function, we can, upon a sequence of rescaling transformations of the Hamiltonian, 
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QP i-^ n QP, obtain a formula involving the ^-function. The fact that this infinite alternate sum of 
transformed Hamiltonians is proportional to the ^-function could be viewed as a confirmation of what 
is expected from (|4.4p and (|4.7p . The characteristic polynomial may then be regarded as a WKB-like 
wave function of the appropriate Hermitian operator. 

In other words, by comparing (|4.12[) to (|4.1ip . we would like to interpret the left hand side as a 
particular limit of the characteristic polynomial det(2: • I — H). In this expression, the Hamiltonian 
is H — ^{QP + PQ), with power-like wave functions ipEiq) — A ^-1/2+*^^^ exact real eigenvalues 
and scaling dimensions ^ — lE. This particular limit of the characteristic polynomial of H is given 
as the rescaling QP n QP, for an integer n, where we introduce the alternate summation over n. 
The claim we make is that in this particular limit the characteristic polynomial of H — \{QP + PQ) 
has zeros at the critical zeros of S,{z). This would be consistent with what we have written in ()4.4p . 

Connes [11] has examined the action of the operator e"^'^, which can be viewed as the generator of 
classical symplectic transformations, albeit on an abstract non-commutative space. In some sense, by 

mapping QP i— > n QP and resumming alternately, we have "fermionized" the generator of the clas- 

00 

sical symplectic transformations. Note also that (— l)"e~"^ can be viewed as a twisted partition 

n=0 

function. As such it can be understood in terms of determinants, as familiar from the physics reinter- 
pretation of index formulas. It is noteworthy that the alternate sums of traces feature prominently in 
Connes's work as well. It would be interesting to understand this relation more precisely. 

This peculiar limit of the characteristic polynomial of a Hermitian operator can be understood 
in direct analogy with what happens in the semiclassical WKB ft — > limit of a simple harmonic 
oscillator. In this case, the exact wave functions (the Hermite polynomials) reduce to the Airy function, 
which can be understood as a WKB limit of the characteristic polynomial of the Hamiltonian for the 
simple harmonic oscillator. The equidistant spectrum of the harmonic oscillator gets rearranged to 
give the zeros of the Airy function in the h ^ Q limit. Note that the zeros of the Airy function have to 
stay coUinear given their origin in the WKB limit of a Hermitian operator. Similarly the continuum 
spectrum of the Hermitian \{QP + PQ) Hamiltonian gets rearranged into the collinear spectrum of 
the zeta function in the limit described in the preceding paragraph and in (|4.12p , which is different in 
detail, but similar in spirit to the canonical WKB limit. 

Finally, we note that the continuous scaling dimension associated with the power law function 
that generates the Mellin transform should be preserved in this limit. This could be responsible for 
the chaotic dynamical properties of critical Riemann zeros, as observed in manifold numerical studies. 
The Fermi-Dirac distribution in (|4.12p might be the evidence for the repulsive ( "fermionic" ) nature 
of the Riemann zeros, which would fit the matrix model description of their correlation properties. 
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